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Abstract 

The infrared limit of the lattice Landau gauge gluon propagator is studied. We 
show that the lattice data is compatible with the pure power law {q^Y*^ solution 
of the Dyson-Schwinger equations. Using various lattice volumes, the infinite 
volume limit for the exponent K is measured. Although, the results allow K = 
0.498 — 0.525, the lattice data favours k ~ 0.52, which would imply a vanishing 
zero momentum gluon propagator. 



1 Introduction and motivation 

The gluon propagator is a fundamental Green's function of Quantum Chromodynamics 
(QCD). Certainly, a good description of this two point function is required to under- 
stand the non-perturbative regime of the theory. Moreover, in the Landau gauge, the 
gluon and the ghost propagators at zero momentum are connected with possible con- 
finement mechanisms. In particular, the Kugo-Ojima confinement mechanism requires 
an infinite zero momentum ghost propagator (2 EJ and the Zwanziger horizon con- 
dition J3l |4] |5l implies either reflection positivity violation or a nuU zero momentum 
gluon propagator From this last condition it can be proved that the zero momentum 
ghost propagator should diverge |6|. 

The computation of the gluon and ghost propagators for the fuU range of momenta 
cannot be done in perturbation theory. Presently, we have to rely either on Dyson- 
Schwinger equations (DSE) or on the lattice formulation of QCD; both are first prin- 
ciples approaches. In the first technique a truncation of an infinite tower of equations, 
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together with the parametrizations of a number of vertices, namely the three gluon and 
the gluon-ghost vertices, is required to solve the equations. In the second, one has to 
care about finite volume and finite lattice spacing effects. Hopefully, the two solutions 
will be able to reproduce the same solution. In this article we are concerned with the 
lattice Landau gauge gluon propagator in momentum space. 



for pure gauge theory. The form factor Z(^^) is the gluon dressing function. 

Recently, in 1 7 1 the DSE were solved for the gluon and ghost propagators in the 
deep infrared region. The solution assumes ghost dominance and is a pure power law 
ISl . with the exponents of the two propagators related by a single number K. For the 
gluon dressing function, the solution reads 



The DSE equations allow a determination of the exponent K = 0.595, which, for the 
zero momentum, implies a null (infinite) gluon (ghost) propagator, in agreement with 
the criteria described above, and a finite strong coupling constant defined from the 
gluon and ghost dressing functions II II II2I . Renormalization group analysis based 
on the flow equation 1131 [141 II5I were able to predict, for the solution (|2ji, the range 
of possible values for the exponent. The results of this analysis being 0.52 < fc < 
0.595, implying, again, a null (infinite) zero momentum gluon (ghost) propagator. A 
similar analysis of the DSE but using time-independent stochastic quantisation II6III7I 
predicted the same behaviour and K — 0.52145. 

The reader should be aware that not all solutions of the DSE predict a vanish- 
ing zero momentum gluon propagator In II8I II9I . the authors computed a solution 
of the DSE with a finite zero momentum gluon propagator Moreover, in 1201 12II it 
was claimed that the gluon and ghost propagators in the deep infrared region are not 
connected via the same exponent K and that there isn't yet a clear theoretical answer 
concerning the value of zero momentum gluon propagator. According to the authors, 
the range of possible values goes from zero to infinity. 

On the lattice, there are a number of studies concerned with the gluon propagator 
1221 1231 l24l 1251 l26l 1271 . On large lattices, the propagator was investigated in l;28,|. 
Although the lattices used had a limited access to the deep infrared region, the au- 
thors conclude in favour of a finite zero momentum gluon propagator, which would 
imply K = 0.5 for the solution (|2ji. In ||29) , the infinite volume and continuum lim- 
its of the Landau gauge gluon propagator were investigated using various lattices and, 
again, the data supported a finite zero momentum gluon propagator. In 1241 . the three- 
dimensional lattice SU(2) Landau gauge propagator was studied with the authors mea- 
suring an infrared exponent compatible with the corresponding DSE solution. 

In order to access the deep infrared region, in we have computed 

the gluon propagator with large asymmetric lattices. The fits to (|2j produced always 
K < 0.5. However, the inclusion of corrections to the leading behaviour given by (|3 
or when the gluon dressing function was modelled, either in the momenta region ^ < 1 




(1) 




(2) 
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GeV or for all the momenta, K becomes larger than 0.5. Moreover, in the prelimi- 
nary finite volume study 1321 1341 of the two point gluon function it came out that K 
increases with the lattice volume and, in that sense, previously computed values should 
be regarded as lower bounds. The previous studies do not provide a clear answer about 
the zero momentum gluon propagator. In this paper, we update the results of our pre- 
vious publications for the infrared region, using a larger set of lattices, that allows a 
better control of the infinite volume extrapolation. Two different extrapolations to the 
infinite volume are performed, namely the extrapolation of K and the extrapolation of 
the propagator, giving essentially the same result. The main conclusions from this in- 
vestigation being that the lattice data is compatible with the DSE solution (|2ji for the 
infrared gluon dressing function, for momenta ~ 150 MeV or lower, and that the infi- 
nite volume extrapolation of the lattice gluon propagator seems to favour a vanishing 
zero momentum gluon propagator 



2 Field Definitions and Notation 

In the lattice formulation of QCD, the gluon fields are replaced by the links 



(3) 



where are unit vectors along direction. QCD is a gauge theory, therefore the fields 
related by gauge transformations 



y{x) U^ix) g\x + ae^), 



(4) 



are physically equivalent. The set of links related by gauge transformations to (x) is 
the gauge orbit of (x). 

The gluon field associated to a gauge configuration is given by 



A^{x + ae^/2) ^ ^ [f/^ (x) - t/^ (x) 



1 



-Tr 



f/^(x)-f/;(x) 



(5) 



up to corrections of order a^. 

On the lattice, due to the periodic boundary conditions, the discrete momenta avail- 
able are 



aLn 



0, 1,...L^-1, 



where is the lattice length over direction pL. The momentum space link is 

X 

and the momentum space gluon field 

X 



(6) 



(7) 



2/^0 



U^,{q)-Ul{-q) 



1 

3^^ 



U^,{q)-Ul{-q) 



(8) 
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The gluon propagator is the gluon two point correlation function. The dimension- 
less lattice two point function is 

{A",iq) Atiq')) = DfAq) V 5{q + q') . (9) 
On the continuum, the momentum space propagator in the Landau gauge is given by 

DfM - S'"^ - ^) D{q') . (10) 

Assuming that the deviations from the continuum are negligable, the lattice scalar func- 
tion D{q^) can be computed directly from ilOi as follows 

^^^'^ = (jV2_iXjV,_l)y I<Tr[A^(^)A^(-g)] ), q^Q, (11) 

and 

= (jV2_^l)jV,y I<T''[-^^(^)^A:(-g)] ): ? = 0, (12) 

where 

Nc = 3 is the dimension of the group, Nd — 4 the number of spacetime dimensions and 
V is the lattice volume. 



3 The Landau Gauge 

On the continuum, the Landau gauge is defined by 

df,A^ - 0. (14) 

This condition defines the hyperplane of transverse configurations 

r = {A: d-A = 0} . (15) 

It is well known [351 that F includes more than one configuration from each gauge 
orbit. In order to try to solve the problem of the nonperturbative gauge fixing, Gribov 
suggested the use of additional conditions, namely the restriction of physical configu- 
rational space to the region 

Q. = {A: d-A = Q, M[A] > 0} C L, (16) 

where M[A] = — V • D[A] is the Faddeev-Popov operator. However, Q. is not free of 
Gribov copies and does not provide a proper definition of physical configurations. 

A suitable definition of the physical configurational space is given by the funda- 
mental modular region A C il, the set of the absolute minima of the functional 




4 



In this article, the problem of Gribov copies' on the computation of the gluon propa- 
gator will not be discussed. For a numerical study on Gribov copies see 1361 1371 and 
references therein. 

On the lattice, the situation is similar to the continuum theory. The Landau gauge 
is defined by maximising the functional 



Fu[g\ = Cf^Re{Tr[g(x)f/^(x)gt(x + A)]} 



x,ii 



where 



Cf = 



1 



(18) 



(19) 



NdNcV 

is a normaUzation constant. Let C/^ be the configuration that maximises Fu[g\ on a 
given gauge orbit. For configurations near on its gauge orbit, we have 

Ft/[l+/a)(x)]=Ft/[l] + ^^/a)"WTr[ X" {U ^{x) - U ^{x ~ fi)) - 



X"{ul[x)-Ul[x^fL)) 



, (20) 



where X" are the Gell-Mann matrices. By definition, f/^ is a stationary point of F, 
therefore 



dF 



iCf 



dco"(x) 



£Tr X"{U^ix)-U^{x~^)) 



X"{u;ix)~u;ix~fi)) 
In terms of the gluon field, this condition reads 



= 0. 



^Tr X"{A^{x + afi/2)-A^{x-afi/2)) 



(21) 



(22) 



or 

£^^A^(^) + (?{a) = 0, (23) 

i.e. i2l\ is the lattice equivalent of the continuum Landau gauge condition. The lattice 
Faddeev-Popov operator M(L') is given by the second derivative of ( I18t . 

Similarly to the continuum theory, on the lattice one defines the region of stationary 
points of ( I18> 

r = {{/ : d-A{U) = 0}, (24) 

'For the 16"' x 128 lattice, the gluon propagator was computed using both the gauge fixing method de- 
scribed here and the gauge fixing method described in 1381 . which aims to find the absolute maxima of Fu [g] 
(see equation 1181 ). For a similar number of configurations, i.e. 164 configurations at j3 = 6.0, we found 
no clear differences in the propagators 1311 . Therefore, in this study it will be assumed that Gribov copies 
do not play a significant role. Note that in 1391 . it was shown that, in the continuum, the expectation values 
measured in CI are free of Gribov copies effects. 
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Lattice 


Update 


therm. 


Sep. 


Conf 


8^ X 256 


70VR+4HB 


1500 


1000 


80 


10^ X 256 


70VR+4HB 


1500 


1000 


80 


12^ X 256 


70VR+4HB 


1500 


1000 


80 


14^ X 256 


70VR+4HB 


3000 


1000 


80 


16^ X 256 


70VR+4HB 


3000 


1500 


155 


18^ X 256 


70VR+4HB 


3000 


1500 


40 



Table 1: Lattice setup used in the study of the gluon propagator volume dependence. 



the Gribov's region Q. of the maxima of 

Q. = {U : d-A{U)=^OmdM{U)>0} (25) 

and the fundamental modular region A defined as the set of the absolute maxima of 
J18> . A proper definition of the lattice Landau gauge chooses from each gauge orbit, 
the configuration belonging to the interior of A. 

In this work, the gauge fixing algorithm used is a Fourier accelerated steepest de- 
scent method (SD) as defined in 1401 . In each iteration, the algorithm chooses 



g{x) = exp 

where 



J %^ ^ (lA_v [C/v W - Ulix)] - trace) 



(26) 



A_v (f/ii = (-^ - «e"v) - f/^i (x) , (27) 

are the eigenvalues of (— <3^), a is the lattice spacing and F represents a fast Fourier 
transform (FFT). For numerical purposes, it is enough to expand to first order the ex- 
ponential in ( I26> . followed by a reunitarization of g(x). On the gauge fixing process, 
the quality of the gauge fixing is measured by 

^ = 4w lTr[A(^)At(x)] (28) 

V JVC ^ 

where 

= L \Pvix-aey) - U\{x) - h.c. - trace] (29) 

V 

is the lattice version of d^xA^ = 0. 



4 Lattice setup 

The gluon propagator was computed using the SU(3) pure gauge, Wilson action, j3 = 
6.0 configurations for the lattices reported in table[l] All configurations were generated 
with the MILC code |http : / /physics . Indiana .edu/^sq/milc. htmll The 
table describes the combined overrelaxed-nheat bath Monte Carlo sweeps, the number 
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of thermalisation sweeps of the combined overrelaxed+heat bath steps, the number of 
combined sweeps separating each configuration and the total number of configurations 
for each lattice. The statistical errors on the propagators were computed using the 
jackknife method. 

The propagators computed in this study have finite volume effects 1321 . An exam- 
ple is seen in figure^where the bare gluon propagator, D{q^) = Z{q^)/q^, is plotted 
for temporal momenta qt and spatial momenta q^ for the smallest and largest lattices. 
Figure|3]shows D{q^) for temporal q <l GeV for the volumes considered in this study, 
including the two lattices 16^ x 48 and 32^ x 64, /3 = 6.0, temporal momenta propa- 
gators from 1281 . For the lattices discussed here, it was observed that the difference 
between the propagators computed using only pure temporal or pure spatial momenta 
becames smaller as the spatial extension of the lattice increases. Moreover, the dif- 
ferences vanish for sufficient high momenta and become larger for smaller momenta. 
Furthermore, it was observed that, for the smallest momenta, the propagator decreases 
as the lattice volume increases, in agreement with what was observed in previous stud- 
ies of the volume dependence 1291 1321 . In this work, we will not discuss what should 
be the right choice of momenta to minimize finite volume effects. Since our smallest 
momenta are pure temporal momenta, the results of the various lattices will be com- 
bined to investigate the volume dependence of the infrared gluon propagator. In the 
following and in order to access the infrared region, only the pure temporal momenta 
will be considered. 



5 The Infrared Gluon Dressing Function 

An analytical solution of the DSE for the gluon dressing function is given by (|2ji. For 
/3 = 6.0, the lattice spacing is a^' = 1.94(5) GeV 1331 . therefore our smallest nonzero 
momentum is 47.6 ± 1 .2 MeV. Since we do not know at what momenta the above solu- 
tion sets in, besides the pure power law, we will consider also polynomial corrections 
to the leading behaviour, i.e. 



with n = 1, 2. In order to be as close as possible to the infrared region, in all fits only 
the smallest range of momenta will be considered, i.e. all reported fits have one degree 
of freedom. The K and /d.o.f. for the various fits are given in table|2] In general, 
the data is well described by any of the above functions. The exception being the 
smallest lattice for a pure power law. K shows finite volume effects, becoming larger 
for larger volumes. Note that for the largest lattices, the various K computed with the 
corrections to the pure power law agree within one standard deviation. Moreover, for 
the largest lattice, the various K agree with each other within less than 1 .5 standard 




(30) 
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Lattice 




Ik 






(H 


-V) 






8^ 


x256 


0.449611 


2 


14 


0.4773 


f37 
-52 


0.02 


0.48271^^ 


0.00 


10^ 


x256 


0.46501^^ 





10 


0.4827 


f49 
-68 


0.25 


0.4765li«4 


0.14 


123 


x256 


0.46631^0 


1 


19 


0.4822 


f51 
-69 


0.21 


0.48491^^ 


0.18 


143 


x256 


0.49181^^ 





09 


0.5053 


f52 
-67 


0.16 


0.49921J™ 


0.06 


16^ 


x256 


0.4859124 





40 


0.5070 


f36 
-50 


0.44 


0.51311^;^ 


1.03 


18^ 


x256 


0.50171^9 





20 


0.5169 


f89 
-70 


0.00 


0.514l}^ 


0.00 



Table 2: K and jd.o.f . from fitting the different lattices. In all fits, the range of mo- 
menta considered starts with the lowest non-vanishing momentum. The errors shown 
are statistical and were computed using the bootstrap method, with the number of boot- 
strap samples being about ten times the number of configurations. For the momenta 
range considered in the fits, the correction associated with a(f- and aq^ + bcf to the 
pure power law are below 20% and 30%, respectively. For the smallest lattice, the 
corrections are larger. However, this lattice is never used in the extrapolations. 



deviations. Figure |3] shows the results of the fits to if as function of the inverse of the 
lattice volume. 

The infinite volume fc can be estimated combining the results of table|2] Assuming, 
for each of the reported fits to the dressing function, either a linear or a quadratic depen- 
dence on 1 /y for Kiy ) and excluding the smallest lattice, it comes out that the figures 
on the first column, the pure power law, are not described by any of these functional 
forms^. Figure |3] includes, besides the k values for each volume, the extrapolated if 
from fitting k:(1 /V) to a linear or quadratic function of 1 jV . Note that the points 
are computed independently for each of the corrections to the pure power law consid- 
ered. The extrapolated K agree within one standard deviation, with the quadratic fits to 
k(\ /y) giving larger errors and lying above the linear fits. For the linear extrapolation 
the results give 

K-oo = 0.5167(55), K-oo = 0.5218(80) (31) 

when the extrapolation uses the data from the quadratic or the quartic corrections to (|2j. 
The corresponding jd.o.f . being 1 .44 and 0.49. For the quadratic fit the extrapolated 
values are 

K-oo = 0.537(13), K-oc = 0.541(20); (32) 

the /d.o.f. being 0.72 and 0.10, respectively. The errors K were computed assum- 
ing gaussian error propagation. The results of the extrapolation to the infinite lattice 

^The smallest y;} jd.o.f . being larger than 5. This applies even if one considers a cubic fitting function. 
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p 




Set I 






Set II 






Set III 






83-18^ 


x256 




10^ 


- 18^ X 256 


123- 


-18^x256 




Lin 


Quad 


Cub 


Quart 


Lin 


Quad 


Cub 


Lin 


Quad 





2.84 


2.60 


1.54 


0.21 


3.08 


1.94 


0.13 


2.99 


0.01 


1 


6.15 


1.95 


2.60 


4.53 


2.66 


2.67 


4.71 


2.92 


5.07 


2 


2.02 


0.44 


0.56 


0.71 


0.68 


0.59 


0.77 


0.72 


0.90 


3 


1.66 


1.02 


1.32 


0.06 


0.85 


1.20 


0.03 


1.04 


0.00 



Table 3: The jd.o.f . for the extrapolations to the infinite volume of the four lowest 
momenta gluon propagator assuming a linear (Lin), quadratic (Quad), cubic (Cub) and 
quartic (Quart) polinomial dependence on 1 jV and for various sets of data. 



volume, point to a value of K in the range 0.51 — 0.56. 

So far we have extrapolated the fc values. Alternatively, one can extrapolate directly 
the gluon propagator and fit the associated dressing function. In figures|3and|5lthe bare 
propagator is plotted, as function of the volume, for the lowest four momenta. For all 
these momenta, we tried linear, quadratic, cubic and quartic extrapolations. In table|3] 
the jd.o.f. is reported for the various extrapolations considered and for the lowest 
four momenta. 

In the following, we will not consider extrapolations using all lattices (set I), be- 
cause it involves our smallest lattice, which is too short ^ 0.8 fm in the spatial direc- 
tions. Note, however, that the data reported in table|3lpoint towards a smooth approach 
to the infinite volume limit, starting from lattices sizes as small as 0.8 fm. Furthermore, 
figures|3and|5lsuggest that we should not use a linear extrapolation. According to the 
figures, a quadratic extrapolation, or a higher power of the inverse volume, seems to 
be more suitable. However, for set II, and for the first non-vanishing momentum, the 
cubic extrapolation gives a quite poor description of the data, when compared to the 
quadratic fit. So, for the reasons explained above, from now on we will consider only 
the quadratic extrapolation using set II and III. 

Figures|4]and|5]show the lattice data together with the two quadratic extrapolations 
to the infinite volume. Note that in both extrapolations, the worst M-o.f. is asso- 
ciated with the first nonzero momentum. Indeed, as can be seen in figure |4] for this 
momentum the data seems to fluctuate more than for the other momenta. For larger 
momenta the behaviour of the gluon propagator as a function of the volume becomes 
smoother. 

In 1291 the infinite volume of the renormalized zero momentum gluon propagator 
was computed. The chosen renormalization condition was 

dr{i\2.,2 = ^; (33) 

the lattice data was renormalized at = 4 GeV. The quoted extrapolated renormalized 
zero momentum gluon propagator being 7.95(13) GeV^^. Our largest momentum 
available is ^ = 3.88 GeV, a slightly lower value. Performing the renormalization in 
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the same way but at the scale fi = 3.88 GeV, it comes 



{ 



6.3 ±1.4 GeV 
10.9 ±0.8 GeV 



2 



set III data, 
set II data. 



Dr{0) 



2 



(34) 



The quoted errors are pure statistical. The value quoted in [29 1 is in between the above 
figures, with the result from the so called set III being almost compatible within one 
standard deviation with 7.95(13) GeV^^ and the result from set II being compatible 
only at three standard deviations. 

The extrapolated propagators can be seen in figure |6j the errors were computed 
assuming gaussian error propagation. For infrared momenta, the extrapolated propaga- 
tors interpolate between the two lattices simulated in f28V Moreover, for zero momen- 
tum the extrapolated propagators are smaller than the propagator of Leinweber et al 
computed with their largest lattice. The extrapolated propagators from fitting the two 
different sets of lattices are compatible, at least, at the level of two standard deviations. 
The propagator computed from set III shows much larger statistical errors. The differ- 
ences between the two propagators are larger for smaller momenta. This difference can 
be used to estimate systematic errors. 

The fit of the dressing function computed from the extrapolated propagator, ob- 
tained from set II, to the first three nonzero momenta to the pure power law gives 
K = 0.5215(29) and has a /d.o.f. = 0.02. Similar fits but using the corrections to 
(|3 produce larger values for the same exponent (~ 0.55). All the computed K"'s are 
compatible with a vanishing zero momentum gluon propagator Moreover, the fit to 
the pure power law produces a number which agrees well with the estimation from an 
extrapolation on k: - see equations (13 U and \32\ . If, instead, we perform the same anal- 
ysis but for the extrapolated propagator obtained from set III, then K = 0.4979(66) and 
/d.o.f. = 0.21 . Again, the fits with corrections to the pure power law give large val- 
ues for K = 0.52 — 0.53. Note that now the K fitting the pure power law is smaller than 
the corresponding value obtained from our largest lattice (see table|2j and is compatible 
with K = 0.5215(29) within two and a half standard deviations. 

6 Results and Conclusion 

We have computed the Landau gauge gluon propagator for large asymmetric lattices 
with different volumes. For each volume, we have checked that the lattice dressing 
function is well described by the DSE solution (|3 for momenta below ~ 150 MeV by 
fitting the above functional form. The exception being our smallest lattice 8^ x 256. 
Note that in all these fits Z{q^)\y=o was never used. In what concerns the infrared 
exponent K, the data shows finite volume effects, with K becoming larger for larger 
volumes. In this sense, all values in table|2can be read as lower bounds on the infinite 
volume figure. 

The infinite volume k was estimated in two different ways: i) from an extrapola- 
tion of the values reported in table |2j assuming a linear and a quadratic dependence 
on the inverse volume; ii) extrapolating the propagator, assuming a quadratic depen- 
dence on 1 /V. The results of the first method are ( 13 H . ( I32> . giving an weighted mean 
value of K" = 0.5246(46). The value for the exponent from extrapolating directly the 
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gluon propagator being K = 0.5215(29), if one considers the five largest lattices (set II), 
and K = 0.4979(66), if one considers the four largest lattices (set III). The first value 
is on the top of the value obtained from extrapolating directly k as function of the 
volume. Moreover, these two values agree well with the time-independent stochastic 
quantisation prediction K = 0.52145 I16II17I and are within the range of possible val- 
ues estimated with renormalization group arguments I13II14|[T51 . 0.52 < K < 0.595. 
All these results point towards a vanishing zero momentum gluon propagator. The 
K = 0.4979(66), obtained from fitting the gluon dressing function, computed from the 
quadratic extrapolation of the gluon propagator to the infinite volume, using only the 
four largest lattices (set III), is compatible with both an infinite and a null zero mo- 
mentum gluon propagator This value is smaller than the measured exponent from 
our largest lattice, K = 0.5017(49), and, although, these two numbers being compat- 
ible within errors, the extrapolated K does not follow the observed behaviour that K 
increases with the lattice volume. This is probably due to the large statistical errors 
observed in the extrapolation using the smaller set of lattices. Finally, one can claim a 
K = 0.498 — 0.525 with the lattice data favouring the right hand side of the interval. 

Although the lattice data seems to favour a null zero momentum gluon propagator, 
the measured extrapolated zero momentum propagator clearly does not vanish - see 
equation ( I34t . As a function of the lattice volume, the bare zero momentum propagator 
becomes smaller for larger volumes. Moreover, in 1321 it was shown that the extrapola- 
tion of the zero momentum propagator, depending on the functional form considered, 
is compatible with both a vanishing and a nonvanishing value. In this study we avoided 
the question of the zero momentum value by fitting only the smallest nonzero mo- 
menta. Solutions to this puzzle require simulations on larger lattices and/or having a 
better theoretical control on the volume extrapolation. Currently, we are engaged in 
improving the statistics for our larger lattices and simulating bigger volumes, aiming 
to improve the infinite volume extrapolations. 
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Figure 1: Bare gluon dressing function for all momenta for the largest (top) and the 
two smallest (bottom) lattices. The figure shows the gluon propagator computed with 
pure temporal momenta qt and pure spatial momenta q^. For the spatial momenta a Z3 
average was performed. 
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Figure 2: Bare gluon dressing function for temporal momenta for all lattices volumes. 
For comparisation, the figure includes the 16^^ x 48 and 32^ x 64 propagators computed 
at 6 =6.0 in l28l. 



15 



o 



2 2k 

(q ) 

(q ) ( 1 + a q ) 

, 2,2k, , 2,4, 
(q ) ( 1 + a q + b q ) 

Linear Extrapolation 

Quadratic Extrapolation 



I 



I 



i 



i 



I 



le-06 2e-06 3e-06 4e-06 5e-06 6e-06 7e-06 8e-06 

i/y 



Figure 3: Jf as a function of the inverse volume. Note that in order to make the plot 
more clear, the K from the fits to the corrections to the pure power law were shifted in 
the 1 /V axis. 
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Figure 4: Bare gluon propagator for the two smallest momenta as a function of the 
inverse volume. For comparison, for zero momentum we include D{0) from the largest 
lattice of EH, i.e. 32^ x 64. 
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Figure 5: Bare gluon propagator as a function of the inverse volume for the third and 
fourth smallest momenta. 
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Figure 6: The bare extrapolated gluon propagator. The full range of temporal momenta 
is seen on the top with a zoom for the infrared region on the bottom. The errors on the 
propagator were computed assuming a gaussian error propagation. For comparisation, 
the figure includes the 16^ x 48 and 32^ x 64 propagators computed at j3 = 6.0 in 1281 . 
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